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Abstract 

Background: Pentavalent antimonials are the first line drugs for the treatment of 
leishmaniasis. Unresponsiveness of Eeishmania spp. to antimonial drugs is a serious 
problem in some endemic areas. Investigations on molecular mechanisms involved 
in drug resistance are essential for monitoring and managing of the disease. Cal- 
cineurin is an essential protein phosphatase for number of signal transduction 
pathways in eukaryotic cells and it has a mediated role in apoptosis. This study 
aimed to determine of biomarker(s) in Glucantime® resiatance strain of E. infan- 
tum. 

Methods: We used cDNA amplified fragment length polymorphism (cDNA- 
AFLP) and real time-RT PCR assays to compare gene expression profiles at the 
mRNA levels in resistant and susceptible E. infantum field isolates. 
Results: The cDNA-AFLP results showed downlegulation of calcineurin in resis- 
tant isolate in comparison with susceptible one. Significant downregulation of cal- 
cineurin (0.42 fold) (P<0.05) was found in resistant isolate compared to susceptible 
one by Real time-RT PCR. 

Conclusion: This is the first report of calcineurin implication in Glucantime® drug 
resistance of field (natural) isolate of E. infantum. Downregulation of calcineurin 
could protect parasites from antimonial-induced apoptosis. 
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Introduction 

Several species of Leisbmama are causative 
agents of a wide spectrum leishmaniasis 
which ranging from self-healing skin 
lesions to dangerous visceral forms (1). Global 
prevalence of leishmaniasis accounts for 12-15 
million people and incidence rate of the disease 
is estimated 2 million new cases annually (2). 
Despite importance of the disease and serious 
attempts of researchers, no desirable vaccine is 
yet available against the disease. Rapid detection 
and appropriate treatment in most clinical forms 
of leishmania-sis have essential role in control of 
the disease (3-5). 

Pentavalent antimonials sbV, such as Glu- 
cantime® have been used clinically as the first 
line drug against leishmaniasis since seven 
decades (6). Although these drugs have been 
administrated long time worldwide, some bio- 
chemical aspects of sbV metabolism in Leish- 
mania are not uncovered (3). The sbV is a 
prodrug and kiUs "Leishmania by sblll which is 
an activated form of sbV. The exact site of 
this activation and involved mechanisms in 
this process are still unknown (7). 

Unfortunately, in recent years the efficacy of 
antimony therapy has been challenged by oc- 
currence of drug resistance. The major prob- 
lem in treatment of leishmaniasis in some re- 
gions such as India is emerging of sbV-resis- 
tant parasites (8). Unresponsiveness to Gluca- 
ntime® has been also reported from Iran (9). 

Various mechanisms are involved in drug re- 
sistance of l^ishmania such as drug entry, drug 
metabolism, drug transport, programmed cell 
death (10). For instance, Aquaglyceroporine 
(AQPl) is a plasma membrane protein of 
l^ishmania and involves in entry of activated 
form of drug sblll to parasite (11). Transfec- 
tion and knock out experiments of AQPl 
have shown that decrease of AQPl led to re- 
sistance (11, 12). TDRl and lmACR2 of para- 
site have been detected which involved in sbV 
reduction (13, 14). It is well documented that 
MRPA is one of the most important gene in 
drug transport and sequestration (12, 15). In 



addition, some studies illustrated the involve- 
ment of Heat Shock Proteins (HSPs) in re- 
sistance by modulating some phase of apopto- 
sis pathway (10). Furthermore, it has been 
suggested that the expression of wild-type cal- 
cineurin could be a key element in the pro- 
cesses of apoptosis. In the eukaryotic cells 
such as cardiac and nerve cells the biological 
function of this protein is dependent to cal- 
cium/caknoduMn (16, 17). 

In spite of the fact that several techniques 
such as RT-PCR, proteomics and microarray 
are currently employed for identification of 
molecular mechanism that are involved in an- 
timonial resistance of l^ishmania, some aspects 
of this issue require further investigations by 
other approaches (18, 19). cDNA-AFLP as a 
method allows researchers to monitor trans- 
criptional changes associated to development 
and modification in cellular functions. In con- 
trast to some techniques such as GeneChip, 
microarray which are relied on specific probs, 
cDNA-AFLP does not need any prior 
knowledge of gene sequences (20, 21). 

In this study we aimed to utilize cDNA- 
AFLP for identification of potential biomark- 
ers in antimonial drug resistance in l^ishmania. 

Materials and Methods 

Isolation of parasites 

Two natural l^ishmania spp. isolated from 
visceral leishmaniasis patients by bone marrow 
aspirations. One of them was clinically re- 
sistant to Glucantime® but another one was 
sensitive to the drug. Resistant case has been 
treated three periods with systemic administra- 
tion of Glucantime® but did not lead to cure. 
Susceptible case after bone marrow aspiration 
has been treated with Glucantime® and led to 
cure. These isolates were confirmed as JL. in- 
fantum by ITSIPCR-RFLP method (22). Para- 
sites were cultured in RPMI 1640 medium 
(Gibco/BRL) supplemented with 10% fetal 
bovine and incubated at 25°C. 
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Ethical Committee of the School of Public 
Health, Tehran University of Medical Sciences 
reviewed and approved the study protocol. 

Amastigote drug susceptibility assay 

Drug susceptibly of amastigotes of -L. infan- 
tum to Glucantime® was determined by culti- 
vation of the parasites in the J774A.1 mono- 
cyte-macrophage mouse line. Briefly, macro- 
phages (5 X lO'* macrophages/well) were cul- 
tured in RPMI1640 with 10% FBS in eight 
chamber LabTek tissue-culture slides and in- 
cubated at 37 C for 24 h. In order to infect 
macrophages by stationary phase promas- 
tigotes, 5 X lO^promastigotes/weU were added 
to macrophages then incubated at 37 C for 4 h. 
Subsequently; cells were incubated for 72 h 
with serial dilutions of Glucantime®. Pentava- 
lent antimony concentrations for sensitive iso- 
late used 2, 4, 6, 8, 10 and 12 |ag/ml and for 
resistant isolate were 35, 40, 45, 50, 55, 60, 65 
and 70 |^g/ml (the doses were used based on 
our previous screening test). Fresh drug was 
added to slides for an additional 72 h. After 
staining by Giemsa, based on counting the 
amastigotes in 100 randomly chosen macro- 
phages, the 1C50 is defined as the effective 
dose of Glucantime® that decreases the sur- 
vival of -L. infantum by 50%. Experience was 
performed triplicate and 1C50 values were de- 
termined by linear regression analysis. 

RNA extraction 

Total RNA was extracted using Tripure kit 
according to the manufactur's protocol (Roch, 
Mannheim, Germany) with minimal modifica- 
tions. Briefly, 1 xlO** promastigotes were 
packed by centrifugation. Subsequently, the 
pellets were lysed in 1ml Tripure reagent and 
followed by adding 200)^1 choloroform. After 
centrifugation, the aqueous phase was col- 
lected and precipitated by adding of 500 \)1 
isopropanol. Then, pellets were treated with 
75% ethanol and air-dried. Precipitated RNA 
was dissolved in RNase free water. The quan- 
tity and qualit}r of RNA was examined using 



nanodrop (ND-1000, Thermo-Scientific Fish- 
er, US) and gel electrophoresis, respectively. 
The extracted RNA was treated with DNasel 
to eliminate any genomic contamination (Qi- 
agen, Germany). 

cDNA AFLP 

cDNA AFLP was conducted as described 
previously (23). Briefly, single strand cDNA 
(sscDNA) was synthesized with of total 
RNA, 10 mM of dNTP (Fermentas, Burling- 
ton, Canada), 20pmol/)al OligodT (Fermentas, 
Burlington, Canada), 20pmol/)^l random hex- 
amer (Fermentas, Burlington, Canada) and 
incubated at 65°C for 10 min following by 
addition of 200 U RevertAid premium Re- 
verse Transcriptase (Fermentas, Burlington, 
Canada) 20 U Ribolock RNase inhibitor (Fer- 
mentas, Burlington, Canada) and A\xL of 5X 
reverse transcriptase buffer containing Tris— 
HCl (PH 8.3) (Fermentas, Burlington, Canada) 
then incubated at 25°C for 10 min following 
by 50°C for 60 min. The integrity of cDNA 
was checked with glyceraldehyde 3-phosphate 
dehydrogenase (GAPDH) gene primers as 
housekeeping (Table 1). The PCR condition 
was: (5min at 94°C; 30 cycles of 30s at 94°C; 
30s at 55°C; 45s at 72°C and 7 min at 
72 °C). Double strand cDNA (dscDNA) was 
generated using 1.5 U Ribonuclease H (Roche, 
Mannheim, Germany), 40U DNA Polymerase I 
(Fermentas, Burlington, Canada), 10 mM of 
each dNTP and incubation at 16 °C for 2 h, 
followed by purification with High Pure PCR 
Cleanup Micro Kit (Roche, Mannheim, Ger- 
many). Double restriction digestion was per- 
formed with 5U Mbo I and 5U EcoR I (Fer- 
mentas, Burlington, Canada) on 5 )ag purified 
dscDNA at 37°C for 2 h. 

Digested dscDNA fragments were Ugated by 
AFLP adaptors (Table 1) 8^g AD EcoR 1, AD 
Mbo I and 4 jxg ad EcoR I, ad Mbo I. Ligation 
was conducted in a final volume of 60 j^l; it 
was done in following conditions: 1 min at 
55 °C, decreasing to 10 °C over 45 min 
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Tablel : Primer and adaptors sequences used in cDNA-AFLP and Real-Time RT-PCR 



NAME 




SEQUENCES 




PRIMERS 


AD ECORl 




5'-ACCGACGTCGACTATCCATGAAG -3' 




Adaptor 


AD ECORl 




5'- AATTCTTCATGG -3' 






AD MBOI 




5' -CACTATCCAGACTCTCACCGCA -3' 






AD MBOI 




5'-GATCTGCGGTGA - 3' 






PR ECORl 


5' 


- ACCGACGTCGACTATCCATGAAGAATTC - 


3' 


Pre-ampliflcation 


PRMBOI 




5'- CACTATCCAGACTCTCACCGCAGATC -3 






SI ECORl 


5'- 


ACCGACGTCGACTATCCATGAAGAATTCC 


-3' 


Sensitive 


S2 ECORl 


5'_ 


ACCGACGTCGACTATCCATGAAGAATTCG 

J. XVj J. V V_J X V_I J. Xv_/ X J. X X / X X VJ J. X VJ J. U. XXX V,_^ V_J 


-3' 




S3 ECORl 


5'- 


ACCGACGTCGACTATCCATGAAGAATTCA 

J. V V,_^ V _^ V_7 _L V V _^ V_7 X V_7 _L X X J- X X V_j X X V_7 J_ XJ. XXwJ J- XJ. XX X X 


-3' 




S4 ECORl 


5'- 


ACCGACGTCGACTATCCATGAAGAATTCT 

J- X v_^ v_j J. X v_j X v_j j_ X X J- X X v_> j. x X _i. xj. x v_j j_ xj. XX x x 


-3' 




SI MBOI 




5'- CACTATCCAGACTCTCACCGCAGATC -3' 






S2 MBOI 




5'- CACTATCCAGACTCTCACCGCAGATC -3' 






S3 MBOI 




5'- CACTATCCAGACTCTCACCGCAGATC -3' 






S4 MBOI 




5'- CACTATCCAGACTCTCACCGCAGATC -3' 






G6PDH F 




5'- ATCAACGACGCACTGCTG-3' 




Housekeeping 


G6PDH R 




5'- TTCATCCGCTTCCTTAGGC-3' 






Calcineurin F 




5'-GTTTTCAGTGGACCCAGGAG-3' 




Target Gene 


Calcineurin R 




5'-TGAAACTGTCGTACACCTTGAA-3' 







(i.e. 1°C per 1 min) then 6U T4 DNA Li- 
gase (Roche, Mannheim, Germany) was added 
to reaction and incubated at 4 °C overnight. 
Purification was carried out by High Pure 
PGR Clean up Micro Kit (Roche, Mannheim, 
Germany) then Pre-ampMfication was per- 
formed with pre-amp primers (Tablel) ac- 
cording to this program: 30 cycles of 94 °C, 30 
s, 64 °C, 30 s, and 72 °C, 45 s and final exten- 
sion at72 °C for 7 min. 

Pre-amplified products were subjected to 
PGR using sensitive primers which containing 
adaptor sequences plus one nucleotide 3' (Ta- 
ble 1). Products resulting from sensitive am- 
plification were separated on 10% non-dena- 
turating polyacrylamide gel electrophoresis 
(PAGE) and stained with silver nitrate 

Isolation, cloning and sequencing 

Selected bands were extracted from the 
PAGE. The eluted DNA was re-amplified by 
selective (sensitive) primers (Tablel) in 30 cy- 
cles with the appropriate selective primers 
(Table 1) from resultant profile. PGR products 
were checked on 1.5% agarose gel then cloned 
in to a pGEM-T Vector System I (Promega, 



Fitchburg. USA). Recombinant plasmids were 
sequenced for Identification of cloned tran- 
scribed-derived factors (TDFs), using univer- 
sal primers (T7 promotor and SP6) (Bioneer, 
Seoul, South Korea). Homology searches were 
done in non-redundant nucleic and protein 
databases BLAST. 

(http://www.ncbi.nlm.nih.gov/BLAST/). 

Real Time PCR analysis 

Real-Time (RT-PGR) was conducted to in- 
vestigate differences in expression patterns of 
identified gene(s) between resistant and sus- 
ceptible isolates. Specific primers were de- 
signed by primer 3 version 0.4.0, (http://- 
frodo.wi.mit.edu/) according to identified tar- 
get gene (Tablel). cDNA was generated from 
Ij^g of total RNA using QuantiTect® Reverse 
Transcription(Qiagen, Germany) according to 
manufacturer's instruction. 

This experiment was done in triplicate with 
20)^1 volumes using IQ SYBR green Super mix 
(Takara, Japan), in an RT-PGR machine (Gor- 
bett, RG-3000, Australia) .The PGR condition 
was as follows: activation at 95°G for 30s, am- 
plification at 95°G for 5 s, 60°G for 34 s for 48 
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cycles. Comparing the cycle thresholds (CTs) 
of the identified genes with housekeeping 
gene (GAPDH) the relative value of the ex- 
pression level of each target gene was analyzed 
by using the 2 



quenced. Sequencing results were searched in 
non-redundant nucleic and protein databases. 



"'^^^ method. 



Results 

Susceptibility of Leishmania infantum 
amastigotes to Glucantime® 

The IC50 values of sensitive isolate was 
5.33+1.15 i^g/ml while this value for resis- 
tance isolate was 43.33 ±5.77 )a,g/ml. 

cDNAAFLP 

cDNA AFLP approach was conducted to 
identify differential transcriptomics expression 
in resistant and suseptible isolates to Glucan- 
time®. Total RNA was extracted from sam- 
ples in logarithmic phase (data not shown). 
The single strand cDNA (sscDNA) was syn- 
thesized as mentioned in materials and me- 
thods. Using spesific primers of housekeeping 
gene (GAPDH), the integrity of cDNA was 
cheked and confirmed (Fig. 1). 



250 bp 



150bp 



50 bp 



Fig. 1: Detection of housekeeping gene, GAPDH, 
(165bp) by reverse transcriptase-PCR (RT-PCR) 
on an ethidium bromide stained agarose gel (1.5%). 
M : 50 bp (base pair) molecular weight marker. S : 
Sensetive, R : resistant, NTC: No template control 

The cDNA AFLP products from different 
primer combinations represented on 8% non- 
denaturating PAGE (Fig. 2). 

To detect involved gene(s) in drug resistance 
mechanisms, different bands containing TDFs 
were excised from PAGEs and cloned into 
TA-cloning vector. The recombinant plasmids 
containing cDNA AFLP products were se- 




SSMboI S2MboI SI EcoRI 
S3 EcoRI S-l EcoRI S2 EcoRI 



M S R S R S R 




Fig. 2: Pattern of TDFs which extracted from 
cDNA- AFLP PAGE. Sensitive amplification of 
cDNA AFLP on a PAGE from three different 
primer combinations: S3 Mbol /S3 EcoRI, S2 
Mbol /S4 EcoRI and SI EcoRI /S2 EcoRI. The 
arrows indicate extracted TDFs with cod number: 
1, 2, 3 (related genes mentioned in table 2).M; (50 
bp) molecular weight marker. S:sensitve. R: 
Resistant 

In this study, three genes including calcineurin 
and two unknown genes were identified with 
different mRNA expression levels in resistant L. 
infantum isolate compared to the sensitive isolate 
(Table 2). 

Validation of cDNA AFLP results using 
teal time RT-PCR 

In order to confirm the cDNA AFLP results, 
we explored expression of detected gene ex- 
pressions in samples. Real time RT-PCR dis- 
covered the downregulation of calcineurin gene 
(Fig. 3). Figure 3 shows significant downre- 
gulation of calcineurin (0.42 folds) (P<0.05) in 
resistant isolate compare to susceptible. 
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Table 2: Differentially expressed transcription-derived fragments identified by cDNA-AFLP 



Code 


Length 


Accession No 


Annotation 


E-value 


No 


(bp) 








1 


120 


XM-001468586.1 


Calcineurin 


2e-31 


2 


180 


XM_001470402.1 


hypothetical proteins 


le-12 


3 


260 


XM_001470402.1 


hypothetical proteins 


3e-13 



Calcinurin 



I 1.2 

I ' 
% 

J 0.6 
S 

* 0.4 
I 0.2 

* 0 




s 

■ R 



Fig. 3: Relative expression pattern of calcinurin 
gene by real time RT-PCR which was measured in 
L infantum antimony sensitive compared to resis- 
tance. The expression of GAPDH was used to 
normalize the data. The values are the mean ± SD 
of at least three different experiments (P<0.05) 

Discussion 

Since development of effective vaccines 
against leishmaniasis has not been successful, 
chemotherapy is the only choice to manage of 
the disease (3). Emerging of large-scale in- 
crease in antimony drugs resistance has been 
appeared as a serious problem in some en- 
demic regions. To address this problem, sev- 
eral laboratories have undertaken different 
strategies (18, 19). In the present study, we 
utilized cDNA-AFLP as a valuable method to 
discovery gene(s) potentially involved in drug 
resistance of JL. infantum. We counted about 
1500 TDFs which were differentially ex- 
pressed in the studied samples. 

We found downregulation of calcineurin in 
Glucantim-resistant isolate compared with 
susceptible. Calcineurin is a calcium and cal- 
moduMn-dependent protein phosphatase and 
the only protein phosphatase regulated by 
Ca ^ and involved in diverse cellular activities 
including, cell survival and apoptosis (24-26). 



Several studies reported the roles of this 
protein in adaptations under different envi- 
ronmental factors, in terms of salt levels and 
temperature changes (25, 26). For instance, in 
yeast, adaptation to high salt environment oc- 
curs by induction of Pmr2p expression via 
calcineurin activation (24). Adaptation to envi- 
ronmental temperature changes by calcineurin 
has been well described in Cryptococcus neoformans 
and Arabadopsis thaliana (25, 27, 28). Odom et al. 
showed its critical role in surviving of C. neofor- 
mans in serum (25). Calcineurin by interaction 
with other biomoleculs such as heat shock pro- 
teins provides suitable theromtolorance and 
virulence in 'Leishmania major (29). 

In spite of the fact that calcineurin directly 
or indirectly is implicated in surviving, a line 
of evidence suggested that under different 
condition it could play adverse function. For 
example, calcineurin could trigger apoptosis in 
different organisms by particular concentra- 
tion of cytosoMc ROS (16, 30). Due to high 
similarities between mammalian and leishma- 
nia calcineurin in terms of function and struc- 
ture(16), we suggest that downregulation of 
calcineurin might have a negative effect on 
apoptosis in leishmania spp. Consistent with 
our findings for first time a study by Dhein et 
al. reported that calcineurin mediates apopto- 
sis in lymphocytes (31). The activation of cal- 
cineurin is depended upon cytoplasmic cal- 
cium ion concentrations (16). In this regard, a 
study using cardiac cells revealed that the in- 
creasing of intracellular Ca^^ levels led to cel- 
lular apoptosis by activation of calcineurin and 
some transcriptional factors (32). Furthermore, 
solid evidence supported the implication of 
elevated intracellular Ca^^ concentrations in 
leishmania parasite cell death (33-35). 
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Apoptosis has been also observed and char- 
acterized by externalization of phosphatidyl- 
serine and nuclear DNA fragmentation in 
amastigotes of L. donovani treatment with an- 
timony compounds (33). Antimony drugs in- 
duced the productions of some oxidative 
agents such as nitric oxide (NO) or hydrogen 
peroxide (H2O2) which play a role in leish- 
manicidal effects of these drugs (36). It is 
quite understandable that oxidative stress is 
responsible for elevation of intracellular Ca ^in 
"Leishmania and due to this event, apoptosis 
could be occurred by calcineurin activation. 

Conclusion 

For the first time, we identified calcineurin 
as a drug resistance target gene in _L. infantum. 
We suggest that downregulation of this gene 
might be involved in survival rate of L^ishma- 
nia by inhibiting the apoptosis. However, fur- 
ther studies are needed to define precise bio- 
logical function of calcineurin in the process 
of l^ishmania drug resistance. 
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